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Abstract: Bats living in close contact with people in Rwanda were tested for evidence of infection with viruses 
of zoonotic potential. Mucosal swabs from 503 bats representing 17 species were sampled from 2010 to 2014 
and screened by consensus PCR for 11 viral families. Samples were negative for all viral families except 
coronaviruses, which were detected in 27 bats belonging to eight species. Known coronaviruses detected 
included the betacorona viruses: Kenya bat coronaviruses, Eidolon bat coronavirus, and Bat coronavirus 
HKU9, as well as an alphacoronavirus, Chaerephon Bat coronavirus. Novel coronaviruses included two 
betacorona viruses clustering with SARS-CoV, a 2d coronavirus, and an alphacoronavirus. 
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Introduction 


Bats are natural reservoirs for a number of pathogens of 
public health concern (Plowright et al. 2015; Shi 2013). For 
example, in Southeast Asia, Pteropus fruit bats are the 
natural reservoirs of the zoonotic paramyxoviruses Hendra 
and Nipah (Chua et al. 2000). In 2002-2003, an epidemic 
of severe acute respiratory syndrome (SARS) caused by a 
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novel coronavirus (SARS-CoV) emerged in China (Drosten 
et al. 2003, Ksiazek et al. 2003), and bats were determined 
to be natural reservoirs and the possible source of the virus 
(Ge et al. 2013; Lau et al. 2010). In 2012, a pathogenic 
paramyxovirus, Sosuga virus, which caused severe illness in 
a patient following contact with bats in Uganda, was sub¬ 
sequently detected in Egyptian fruit bats (Rousettus aegyp- 
tiacus) (Amman et al. 2015). In Africa, Egyptian fruit bats 
(R. aegyptiacus) are reservoirs of Marburg virus (Towner 
et al. 2008), and antibodies against Zaire ebola virus have 
been detected in the same species (Pourrut et al. 2009), 
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while antibodies against Bombali ebolavirus have been de¬ 
tected in Little free-tailed bat ( Chaerephon pumilus ), An¬ 
golan free-tailed bat ( Mops condylurus) (Goldstein et al. 
2018). 

To assess the risk presented by human-bat contact, the 
USAID Emerging Pandemic Threats PREDICT project has 
been conducting viral surveillance in wildlife in more than 
35 countries to detect viruses of zoonotic potential, 
including in Rwanda (https://ohi.vetmed.ucdavis.edu/prog 
rams-projects/predict-project). 

Rwanda is continental Africa’s most densely populated 
country (Butler 2004), and bats are frequently observed in/ 
around urban centers and adjacent to wildlife protected 
areas. Ecotourism centered on mountain gorillas in Vol¬ 
canoes National Park also drives tourism to the nearby bat 
roosting “Musanze Caves” (Spenceley et al. 2010; Joachim 
2013). 

Between 2010 and 2014, biological sampling of bats 
was conducted at urban and rural sites in Rwanda char¬ 
acterized by an intense human-wildlife interface, including 
in and around Volcanoes National Park and the Musanze 
Caves. Anthropogenic activities around sampling sites were 
classified according to human livelihoods and activities, 
including ecotourism, crop farming, and national parks. 

A total of 503 bats belonging to 17 species were cap¬ 
tured at 25 sites (Fig. 1) following established PREDICT 
protocols for bat capture and sampling (PREDICT 2017) 
(Table 1). Following capture, bats were photographed, 
measured, and identified to species level as close as possible 
(Kingdon et al. 2013; Patterson and Webala 2012). Date, 
site name, season, apparent species, sex, age class (deter¬ 
mined by the degree of epiphyseal-diaphyseal fusion) 
(Anthony 1988), reproductive status, and mass data were 
recorded. For handling and restraining larger fruit bats, 
light anesthesia was induced using inhalational isoflurane 
and oxygen (Fluriso™, Teva UK, Limited, Castleford, UK). 
Smaller fruit bats and insectivorous bats were physically 
restrained during sampling. All bats were released at the 
capture site within 3 h of capture. 

Sampling for viral family screening included collection 
oral and rectal mucosal or fecal swabs. Swabs were placed 
into viral transport media (BD Universal Viral Transport™ 
medium, Becton, Dickinson and Co., Sparks, Maryland) 
and NucliSENS Lysis Buffer (bioMerieux SA™, Marcy 
l’Etoile, France) and stored in liquid nitrogen for transport 
and transfer to a — 80 °C freezer. 

We extracted RNA from oral and rectal swab samples 
using the Qiamp Viral Mini kit™ (Qiagen Inc., Valencia, 


CA) and reverse transcribed into cDNA using Superscript 
III™ (Invitrogen Corp, Carlsbad, CA). Primers targeting 
the housekeeping gene /Lactin were used to ensure the 
presence of amplifiable nucleic acid in RNA extracts 
(Goldstein et al. 2004). Samples were screened by consensus 
PCR targeting conserved gene fragments using established 
assays know to detect corona (Quan et al. 2010; Watanabe 
et al. 2010), alpha (Sanchez-Seco et al. 2001), arena (Lo¬ 
zano et al. 1997), bunya (Briese et al. 2007), filo (Zhai et al. 
2007), flavi (Moureau et al. 2007), hanta (Raboni et al. 
2005), influenza (Anthony et al. 2012), paramyxo (Tong 
et al. 2008), lenti (Courgnaud et al. 2001), and rhabdo 
(Wray et al. 2016) viruses. Bands of the expected size were 
excised from 1% agarose and purified using the Qiaquick 
kit (Qiagen Inc.). Purified PCR products were cloned 
(pCR4-TOPO vector; Invitrogen Corp.) and sequenced 
(ABI 3730 Capillary Electrophoresis Genetic Analyzer; 
Applied Biosystems, Inc., Foster City, CA). Sequences were 
analyzed and edited using Geneious (Version 6.0.3) and 
compared with known sequences in the Genbank database. 
Species identification of PCR-positive individuals was 
confirmed by DNA bar coding of the cytochrome b ( Cytb) 
and cytochrome oxidase subunit 1 (C07) mitochondrial 
genes (Townzen et al. 2008). Fisher’s exact test was used to 
examine the association of viral positivity with age and 
season using STATA 13.0 software. The level of significance 
was set at P < 0.05 (Raymond and Rousset 1995). 

No alpha, arena, bunya, filo, flavi, hanta, influenza, 
paramyxo, lenti, or rhabdo viruses were detected in oral or 
rectal swabs. Coronaviruses (CoV) were detected in 27 
(5.4%) of the 503 bats sampled. Twenty-two of the 27 
coronavirus positive bats belonged to three species: Straw- 
colored fruit bat ( Eidolon helvum; 13 positives of 111 
sampled, 11.7%), Geoffroy’s horseshoe bat ( Rhinolophus 
clivosus; 5 of 24, 24.2%), and Egyptian fruit bat ( Rousettus 
aegyptiacus; 4 of 36, 11.1%). Coronavirus positive bats were 
sampled at 11 different sites, including in the Musanze 
Caves (Table 1). Subadult bats were more likely to be 
positive for CoV than adults (6/36 vs. 21/429; P = 0.04). 
No coronaviruses were detected in juvenile bats (n = 11). 
There were no observed differences between seasons (Dry/ 
Rainy; 3/90 vs. 24/386; P = 0.445). 

Coronavirus sequences were classified as belonging to 
different viral species according to established cutoffs and 
methods (Anthony et al. 2017b). We detected four known 
coronaviruses (Table 2) and four new coronaviruses (Ta¬ 
ble 3; Fig. 1). One of the new betacoronaviruses, Coron¬ 
avirus PREDICT_CoV-43, was detected in Hipposideros 
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Figure 1. Coronaviruses phylogenetic analysis performed in Anthony et al. (2017b). Sequences are collapsed into clades, representing operating 
taxonomic units (sequences sharing equal or more than 90%). The sequences included in this manuscript are indicated in the boxes. 


ruber and R. clivosus bats co-roosting in bat tourism caves 
(Site 10; Fig. 1; Table 3). Comparison of the conserved 
polymerase gene fragment sequences to other known 
coronaviruses indicated that Coronavirus PREDICT_CoV- 
43 clustered near the SARS-like coronaviruses but suggests 
it may be a distinct virus based on the conserved fragment 
sequence, as it showed only 84% nucleotide similarity to 
SARS-CoV (Genbank accession no. NC_009694). The 
second new betacoronavirus, PREDICT_CoV-44, was de¬ 
tected in two Hipposideros caffer bats trapped in Nyungwe 
National Park (Site No. 11) and in a R. clivosus bat in 
tourism caves at Site No. 10 (Fig. 2). Although the con¬ 
served sequence fragment also clustered with other beta- 
coronaviruses, it was quite divergent, showing only 79% 
nucleotide similarity to others in the group. 

The 2d betacoronavirus, PREDICT_CoV-66, was de¬ 
tected in one Rousettus angolensis bat in Nyungwe National 
Park (Site 11; Fig. 2) and showed 84% nucleotide similarity 


to the closest recognized coronavirus, Kenya bat Coron¬ 
avirus BtKY84 (Genbank accession no. GU65428) found 
previously in E. helvum. The only alphacoronavirus PRE- 
DICT_CoV-42 was detected in a R. clivosus bat in tourism 
caves (Site No. 10). This virus sequence showed only 85% 
nucleotide similarity to the closest recognized coronavirus, 
Kenya bat Coronavirus BtKY69 (Genbank accession no. 
GU65413), found previously in horseshoe bats ( Rhinolo - 
phus species). 

Phylogenetic analyses of complete genome sequences 
of coronaviruses from bats, humans, and other vertebrates 
suggest that bats may be the reservoir hosts from which all 
coronavirus lineages originated (Vijaykrishna et al. 2007; 
Anthony et al. 2017a), and several studies document the 
diversity of bat coronaviruses globally (Dominguez et al. 
2007; Annan et al. 2013; Anthony et al. 2017b). 

In this study, sequences representing two novel coro¬ 
naviruses that clustered with the SARS-like coronaviruses 
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Table 1. 

CoV RNA. 

Bat Species Sampled at 25 Sites 

in Rwanda for 11 Viral Families, with Numbers and Percentages of Bats Testing Positive for 

Site no.* 

Species 

Risk interface 

Season 

Total tested (oral and 
rectal swabs) 

No. of positive 

% Positive 

2 

Epomophorus labiatus 

Homestead 

Rainy 

17 

0 

0 

25 

Mops condylurus 

Homestead 

Rainy 

78 

0 

0 

18 

M. condylurus 

Homestead 

Dry 

20 

0 

0 

12 

Hipposideros caffer 

Ecotourism site 

Rainy 

15 

1 

7 


Rhinolophus clivosus 

Ecotourism site 

Rainy 

19 

4 

21 


Rousettus aegyptiacus 

Ecotourism site 

Rainy 

9 

0 

0 


Otomops martiensseni 

Ecotourism site 

Rainy 

1 

0 

0 

7 

R. clivosus 

Ecotourism site 

Rainy 

5 

1 

20 


R. aegyptiacus 

Ecotourism site 

Rainy 

13 

0 

0 


Hipposideros ruber 

Ecotourism site 

Rainy 

2 

0 

0 


Nycteris hispida 

Ecotourism site 

Rainy 

1 

0 

0 

3 

Eidolon helvum 

Ecotourism site 

Rainy 

53 

4 

7,5 

16 

E. helvum 

Homestead 

Dry 

15 

0 

0 

8 

E. helvum 

Homestead 

Rainy 

9 

2 

22 

5 

E. helvum 

Crop farming 

Rainy 

9 

6 

67 

6 

R. Aegyptiacus 

Ecotourism site 

Rainy 

9 

2 

22 

11 

Neoromicia tenuipinnis 

National park 

Rainy 

5 

0 

0 


Myonicteris angolensis 

National park 

Rainy 

6 

1 

17 


H. caffer 

National park 

Rainy 

1 

0 

0 

23 

Stenonycteris lanosus 

National park 

Rainy 

1 

0 

0 


M. angolensis 

National park 

Rainy 

23 

0 

0 

10 

Epomophorus labiatus 

Crop farming 

Rainy 

6 

0 

0 


S. lanosus 

Crop farming 

Rainy 

5 

0 

0 


R. Aegyptiacus 

Crop farming 

Dry 

5 

2 

40 

4 

E. helvum 

Ecotourism site 

Rainy 

5 

1 

20 

22 

N. hispida 

Fishing area 

Rainy 

1 

0 

0 


E. Labiatus 

Fishing area 

Rainy 

3 

0 

0 


Lavia frons 

Fishing area 

Rainy 

1 

0 

0 

21 

N. hispida 

Ecotourism site 

Rainy 

1 

0 

0 


E. Labiatus 

Ecotourism site 

Rainy 

2 

0 

0 

24 

M. condylurus 

Homestead 

Rainy 

20 

0 

0 

20 

Chaerephon pumilus 

National park 

Rainy 

8 

0 

0 


M. condylurus 

National park 

Rainy 

12 

0 

0 

9 

E. Labiatus 

Homestead 

Rainy 

5 

0 

0 

13 

E. Labiatus 

Ecotourism site 

Rainy 

35 

0 

0 


N. tenuipinnis 

Ecotourism site 

Rainy 

1 

0 

0 


Neoromicia cf zuluensis\ 

Ecotourism site 

Rainy 

2 

0 

0 

17 

E. helvum 

Homestead 

Rainy 

20 

0 

0 

19 

Scotophilus viridis 

Homestead 

Rainy 

4 

0 

0 


C. pumilus 

Homestead 

Rainy 

3 

1 

33 


E. Labiatus 

Homestead 

Rainy 

18 

0 

0 

15 

E. Labiatus 

Homestead 

Rainy 

2 

0 

0 


S. viridis 

Homestead 

Rainy 

1 

0 

0 




Coronaviruses Detected in Bats in Close Contact with Humans in Rwanda 


Table 1. 

continued 






Site no.* 

Species 

Risk interface 

Season 

Total tested (oral and 
rectal swabs) 

No. of positive 

% Positive 

1 

M. Angolensis 

Homestead 

Rainy 

16 

0 

0 


E. Labiatus 

Homestead 

Rainy 

9 

1 

11 


S. viridis 

Homestead 

Rainy 

1 

0 

0 

14 

N. hispida 

Homestead 

Dry 

6 

0 

0 

Total 




503 

27 

54 


*Bats sampling sites can be visualized in Figure 2 by site number. 


Table 2. Known Coronaviruses Detected in 

Bats by Species and Specimen Type. 



Bat species 

c-PCR 

Sample 

Year of collec- 

Site 

Risk 

Virus name 

Genbank 


positive 

tested 

tion/ 

no 

interface 


No. 




season 





Epomophorus labiatus 

1 

Rectal 

2013/Rainy 

1 

Home stead 

Strain of Kenya bat coron- 

KX285830 



swab 




avirus/BtKY56/BtKY55 


Chaerephon pumilus 

1 

Rectal 

2013/Rainy 

6 

Ecotourism 

Strain of Chaerephon 

KX285828 



swab 



site 

bat / coronavirus/Kenya/ 
KY22/2006 


Eidolon helvum 

4 

Rectal 

2012/Rainy 

3 

Ecotourism 

Strain of Eidolon bat coron- 

KX285106 



swab 



site 

avirus/Kenya/KY24/2006 


Eidolon helvum 

1 

Rectal 

2012/Rainy 

4 

Ecotourism 

Strain of Eidolon bat coron- 

KX285107 



swab 



site 

avirus/Kenya/KY24/2006 


Eidolon helvum 

6 

Rectal 

2013/Dry 

5 

Crop farm- 

Strain of Eidolon bat coron- 

KX285108 



swab 



ing 

avirus/Kenya/KY24/2006 


Rousettus aegyptiacus 

2 

Rectal 

2013/Rainy 

6 

Ecotourism 

Strain of Bat coronavirus 

KX286259 



swab 



site 

HKU9 


Rousettus aegyptiacus 

2 

Rectal 

2013/Rainy 

1 

Home stead 

Strain of Kenya bat coron- 

KX285819 



swab 




avirus/BtKY56/BtKY55 


Eidolon helvum 

2 

Rectal 

2012/Rainy 

3 

Ecotourism 

Strain of Eidolon bat coron- 

KX285822 


19 

swab 



site 

avirus/Kenya/KY24/2006 



were detected in bat tourism caves and other sites where 
people and bats come into close contact in Rwanda. One 
virus (PREDICT CoV-43) was detected in both SundevalFs 
roundleaf bat ( H . ruber ) and Geoffrey’s horseshoe bat (R. 
clivosus) that were co-roosting in the Musanze Caves (Site 
No. 10). The high sequence similarity of the viral fragment 
detected in both bat species suggests that this virus may 
have the ability to be maintained in more than one host or 
that cross-species transmission may occur. Studies have 
found that viral sharing and cross-species transmission may 


be important factors that contribute to emergence of novel 
coronaviruses and recombination of bat coronaviruses (Lau 
et al. 2010; Johnson et al. 2015). 

While the known coronaviruses detected in this study 
have been identified in other geographical areas and in 
different bat species (Tao et al. 2012; Drexler et al. 2010), 
we report their first detection in Rwanda. The bat coron- 
avirus HKU9 was previously detected in Rousettus lesche- 
naulti bats in China (Tang et al. 2006), and now a strain of 
this virus has been detected in bats in Rwanda. Similarly, 
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Table 3. Novel Coronaviruses Detected in Bats by Species and Specimen Type. 



Bat species c-PCR Sample tested Year of collection/ Site no. Risk interface 

positive season 

Virus name 

Genbank no. 


Rhinolophus clivosus 

1 

Oral swab 

2011/Rainy 

12 

Ecotourism sites 

PREDICT, 

_CoV-43 

KX285821 

Rousettus angolensis 

1 

Rectal swab 

2013/Rainy 

11 

National park 

PREDICT, 

_CoV-66 

KX285426 

Hipposideros caffer 

1 

Rectal swab 

2011/Rainy 

11 

National park 

PREDICT, 

_CoV-44 

KX285826 

Rhinolophus clivosus 

1 
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Figure 2. Map of Rwanda showing the bats sampling sites for corona viral surveillance during the study. 


Kenya bat coronavirus/BtKY56/BtKY55 in R. aegyptiacus , 
Chaerephon bat coronavirus/Kenya/KY22/2006 in Chaere- 
phon pumilus , and Eidolon bat coronavirus/Kenya/KY24/ 
2006 in Eidolon helvum were first detected in Kenya in 2006 
(Tao et al. 2012). We report the presence of these viruses in 
these same bat species in Rwanda, indicating a wider geo¬ 
graphic distribution of these viruses in Eastern Africa, likely 


due to the widespread distribution of their bat hosts 
(Drexler et al. 2010; Gloza-Rausch et al. 2008). 

In conclusion, bats in Rwanda carry novel and known 
coronaviruses, a family of viruses from which novel viruses 
have caused human pandemics. However, bats play 
important ecological roles and their elimination as a con¬ 
trol measure is not recommended or warranted. We rec- 
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ommend additional surveillance and longitudinal studies to 
further understand the ecology of bat coronaviruses and 
the extent of human-bat interactions to identify strategies 
for public health protection and bat conservation. 
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